Theory of spin-Hall transport of heavy holes in semiconductor quantum wells 
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Based on a proper definition of the spin current, we investigate the spin-Hall effect of heavy 
holes in narrow quantum wells in the presence of Rashba spin-orbit coupling by using a spin-density 
matrix approach. In contrast to previous results obtained on the basis of the conventional definition 
of the spin current, we arrive at the conclusion that an electric-field- induced steady-state spin-Hall 
current does not exist in both, pure and disordered infinite samples. Only an ac field can induce a 
S ' spin-Hall effect in such systems. 
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^j- ' I. INTRODUCTION 

(N; 

The emerging field of spintronics stimulates extensive studies of the spin-Hall effect (SHE) in semiconductors with 
spin-orbit coupling. Recently, an intrinsic SHE, which is entirely due to spin-orbit coupling, has been predicted to 
appear in two-dimensional (2D) electron systems Q and in p-doped bulk semiconductors [J- It has been argued 
,-U ■ that this effect occurs even in the absence of any scattering events and that the spin-Hall conductivity is given by a 
universal constant. Subsequently, a detailed treatment of elastic scatteringrevealed that scattering effects can lead to 
a complete cancellation of the total intrinsic spin-Hall current @j El H III Li • Especially the simple Rashba model of 
a 2D electron gas with its unusual properties gave rise to many controversies. Other models that have been studied 

S seem to be more robust against impurity scattering [8j • This discussion is nowadays replaced by a more fundamental 

debate about the proper definition of the spin current [2, Ei fill . H21 Ha . Il4| . Due to non-conservation of spin, which 
results from spin precession, it has been argued that its definition is largely a matter of convenience |15l | . In analogy 
£h ' to the charge current, most researchers identified the spin current simply with the expectation value of the product 
of spin and velocity observables. Based on this definition, the spin current is expressed from a technical point of view 

i ^**i ■ by particular elements of the spin-density matrix namely fy (k, k | t) \ K= o (with the spin quantum numbers A, A', the 
time variable t, and the quasi-momenta k — (fej + fea)/2, k = k\ — k^), which depend only on the quasi-momentum 

t-H k. These quantities are conventionally calculated from multiband Boltzmann equations. However, this procedure 

is not sufficiently general and even fails for particular charge transport problems. There is a subtlety related to 
this definition of transport coefficients. Already for the charge transport, this procedure is only justified, when the 
interaction Hamiltonian of the system commutes with the position operator. Otherwise one has to go back to the 
more general definition, which expresses the charge current through the time derivative of the dipole operator |16lll7| . 
In this case, the current of the multiband system is not obtained from the special elements of the density matrix 

\Q> f\>(k I t)> but is calculated from quantities V K f\i(k, k | t) | K= o> which requires the consideration of the k dependence 
in kinetic equations. For charge transport, the theory of small polarons provides a famous example. In accordance 
with these experiences in the charge transport theory, a proper definition of the spin current has recently been put 
forward 11]. According to this proper definition, the spin current is expressed by the time derivative of the spin 
displacement. This concept of spin current has a number of remarkable advantages 11]. It is in accordance with the 
near-equilibrium transport theory, satisfies the Onsager relations, and provides vanishing spin currents for Anderson 
insulators. As a result, by determining conjugate forces, thermodynamic and electric measurements of the spin current 
become possible. To further illustrate the advantage of the proper definition, let us consider a particular feature of the 
O ■ spin current. Due to its symmetry with respect to time inversion, a finite spin current can exist even in equilibrium 
Q. For simplicity we treat the Rashba model of a 2D electron gas. According to the conventional definition, there is 

• tH , a field-independent stationary spin-Hall current that is completely independent of the spin accumulation. In contrast, 
the spin-Hall current according to the proper definition turns out to be due to the initial time variation of the spin 
polarization and disappears in the steady state, when the spin polarization becomes constant |l3| . 

In the present work, we study the SHE for heavy holes in III-V semiconductor quantum wells on the basis of the 
proper definition of the spin current by systematically deriving and solving the kinetic equations for the spin-density 
matrix. Recently, experimental observations of the SHE have been reported for such a system [l8| . The analysis 
8] of the experimental results seems to provide evidence for a close correspondence between the detected edge spin 
accumulation and the bulk spin currents as described by the conventional spin Hall theory. However, this conclusion 
cannot be considered to be final as the spin polarization near the edges depends on the boundary conditions and may 
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not be directly induced by spin-Hall currents far from the boundaries. The treatment of inhomogeneous and/or finite 
samples requires the calculation of fy(k, K \ t), the k dependence of which introduces a strong coupling between the 
spin and charge degrees of freedom even in the absence of an electric field. The theory of spin transport in finite 
samples has its particular own challenges. 

II. THEORY 

We focus on a model for the lowest heavy hole subbands, which is described by a one-particle Hamiltonian for heavy 
holes in narrow quantum wells being subject to spin-orbit interaction of the Rashba type that results from structural 
inversion asymmetry. The Hamiltonian of the 2D hole gas in the second quantized form 

H = ^ a l\ \ £ k - £f] afcA - J^ (ftw fc • S\\')a} kX a k x> 

fc,A fc,A,A' 

- ieEj2 V K 4_fA a fc+§A| K=o + u^H a fcA a fc'A, (1) 

k,X k,k f A 

is composed of the creation (a kx ) and annihilation (a k \) operators with quasi-momentum k = (k x , k y , 0) and spin A. 
The dispersion relation of free in-plane motion is given by e k = h 2 k 2 /2m. Throughout, the Fermi energy ef is always 
assumed to be the largest relevant energy scale. The strength of the "white- noise" disorder scattering is given by the 
scattering rate H/t = 2ttupf with pp and u denoting the 2D density of states in the absence of spin-orbit coupling 
and the scattering strength, respectively. For simplicity, the electric field E is oriented along the x axis. The cubic 
Rashba spin-orbit coupling 19, 20, 21] is obtained from the Pauli matrices and the energy 

ha k = ^[ i (kl-k 3 _),(kl + k 3 _),0], (2) 

with k± — k x ± iky and fvjJk = ak 3 . a denotes the spin-orbit coupling constant. 

All information that is necessary to determine the kinetic observables are contained in the spin-density matrix, the 
elements of which are grouped in the following form 

/(fc,«|t)=^/ A A (fc,K|i), /(fc,K|t) = X)/A'(*>«l*)*AA'. (3) 

A A, A' 

Although we treat a homogeneous 2D infinite hole gas, it is indispensable to retain both the k and k dependence of 
the spin-density matrix in order to calculate the proper spin current. In the thermodynamic equilibrium [E = 0) and 
to lowest order in a, the spin degree of freedom is characterized by the vector 

f eq (k | t) = -hui k ^^, (4) 

Oe k 

with n(s k ) being two times the Fermi distribution function. With respect to the time dependence, we prefer the 
consideration of Laplace-transformed kinetic equations, which are derived by exploiting the Born approximation and 
by keeping only the lowest-order contributions of the spin-orbit interaction in the collision integral. Due to the n 
dependence, the resulting kinetic equations 

sf - -(« • k)f + iuj K {k) ■ /+ ^V k f = -(/- /) + n(s k ), (5) 

m h t 

m h 

t t de k t de k de k 

are coupled to each other. An integration over the angle of the vector k is indicated by a bar over the respective 
quantity. In addition, we make use of the abbreviation 



Next, the spin current is treated. According to Refs. [llLIT^| . the proper spin current is calculated from the derivative 
of the spin displacement 
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3i(t) = -^-frYl Yl r m,t(al nXi a m \ 2 )t^\ 1 \ 2 , 



(8) 
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where r m is the position operator and the factor 3S./2 refers to the angular momentum of the heavy holes. Performing 
a Fourier transformation and an integration by parts, the Laplace transformed spin-Hall current is expressed by the 
equivalent equation 



f v {s) = -is— Y, d^f Z ( k > K \ s ) 

k v 



(9) 



K=0 



Within the linear response approach, we need only the first-order corrections of the spin-density matrix with respect 
to the electric field E and the quasi-momentum n. Based on perturbation theory, analytical results for this quantity 
are derived. Details of the calculation are given in the Appendix. 

To start our analysis, let us first treat the field-induced spin accumulation. From the analytical solution in Eq. (A. 4) 
together with (A. 5) and (A. 6), we conclude that there is no field- induced spin accumulation in the cubic Rashba model 

(foE ~ 0). This result should be compared with the finite current-induced spin polarization in the linear Rashba 
model [23|. Furthermore, in contrast to the linear Rashba model there is no spin current, which is independent of the 
electric field and which exists even in thermodynamic equilibrium. Both peculiarities of the cubic Rashba model are 
closely related to each other and are compatible with recent studies based on diffusion equations 23] . 

To calculate the spin-Hall current, we use the analytical results from the Appendix and obtain for its nonzero 
component an expression 
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which indicates that all occupied states below the Fermi energy contribute. The very same feature, which is also 
observed in the linear Rashba model, seems to be a generic property of the intrinsic SHE [131 |21| . In Eq. l|l(Jf) . 
the abbreviation a = s + 1/r is used. Integrating by parts, we obtain another equivalent result, in which at zero 
temperature (T = 0) only states on the Fermi surface play a role 



Jy( s ) = o— a 
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The time dependence of the spin-Hall current is obtained from this equation by applying an inverse Laplace transfor- 
mation 



jy(t) = g — exp (-- ) ^2n'(e k )(uj k t)sm(uj k t). 



(12) 



After the electric field is switched on at time t = 0, jy(t) exhibits damped oscillations. In the steady state, there is 
no spin-Hall current. At zero temperature, our general result takes the simple form 
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2tt h 2 (<j 2 +4a 2 k 6 F ) 2 ' 



(13) 



where hp denotes the Fermi momentum. From this equation, we obtain for the spin-Hall conductivity of a perfect 
crystal (r — > oo) the result 



a S H (w — > 0) 



9e / u> 
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(14) 



which vanishes in the limit u> — > 0. This conclusion completely contradicts previous findings derived on the basis of 
the conventional definition of the spin current |20|. We regain this published result for the spin-Hall current from 
Eqs. (A. 4) to (A. 6). From the general expression for T = 



J» 



9 eEa 2 



fvp 



2tt H 2 s a 2 



4a 2 k e F ' 



(15) 



the universal spin-Hall conductivity 



a sH ( W ^0) = - J£j, (16) 



is obtained. These results, which have previously been derived by an alternative approach 20], strongly deviate from 
Eqs. (|13l) and i|14|) . According to the widespread reasoning, the SHE is present and robust against disorder in the 
cubic Rashba model. In contrast, based on the proper definition of the spin current, we come to the conclusion that 
the SHE is absent in an infinite heavy hole gas. 

III. SUMMARY 

Recently, it has been recognized that the SHE has an intrinsic contribution due to spin-orbit interaction in a perfect 
crystal. This assertion has occupied a great deal of attention because of its potential for electronic devices with low 
power consumption. A theoretical controversy about the disorder effect that can completely eliminate the intrinsic 
SHE seems to be settled now: the intrinsic SHE is absent only in the simple model of a 2D electron gas with Rashba 
spin-orbit coupling that turns out to exhibit anomalous properties. In all other generic systems that have been 
specifically studied, the SHE is present and robust against disorder Q. However, this conclusion is derived on the 
basis of an approach for the SHE that revealed rather unconventional properties |9l Hoi ITU Il3. ITa . ll4 J so that serious 
doubts arose on its physical relevance. A recent proper definition of the spin current |lll Il3l | resolves a number of 
difficulties of former approaches and is in line with Onsager relations that allow the application of the near-equilibrium 
transport theory. On the basis of this physically motivated concept of spin transport, we arrive at the conclusion 
that neither in the linear [Tj| nor in the cubic Rashba model an intrinsic SHE exists. Only an ac electric field is 
able to induce a spin-Hall current. This conclusion does not contradict recent experiments |TSj . which revealed edge 
spin accumulations. The observed spin polarization near the boundaries certainly depends to a large extend on the 
boundary conditions and may not be simply due to spin Hall currents in the bulk. In an inhomogeneous system, there 
is always a strong coupling between the charge and spin degrees of freedom that can give rise to a number of generic 
effects not present in an infinite homogeneous 2D electron or hole gas. It is an interesting future task to extend the 
approach applicable for the bulk to the spin-transport phenomena in finite systems. 
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APPENDIX A: SOLUTION OF KINETIC EQUATIONS 

The proper definition of the spin current in Eq. Q sets up our calculational scheme as an perturbational approach 
with respect to n. In addition, we restrict ourselves to the linear response regime, where only the lowest-order 
contribution in the electric field is taken into account. The formal solution of the kinetic Eq. © has the form 

r_ <?X -2tDfc x X + AuJk(uik • X)/a 
a 2 + AujI 

where a = s + 1/r. X comprises all angle integrated quantities as well as the k and field contributions. This result 
can be used to execute the perturbational approach step by step. First, we treat the equations for k — and E = 
and obtain immediately 

foo = = foo, foo = -fu3h — , foo = 0, (A2) 

s s 

where the indices 00 refer to the order in k and E. n' is a short-hand notation for dn(ek)/dek. Next, the lowest-order 
correction due to the electric field is calculated (E ^ 0, K = 0). The kinetic Eq. (0 for the charge degree of freedom 
is easily solved 

f 0E = - ——n', 7^=0. (A3) 

as m 



Based on Eq. (A.l), we obtain for the spin contribution 

aR aE - 2u k x Rq E +4w fc (u fc • R 0E )/<J 
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(A6) 



From this solution, we conclude that there is no field-induced spin accumulation f 0E = 0. This peculiarity is specific 
for the cubic Rashba model. Within the framework of the conventional theory, all components of the spin-density 
matrix that determine the spin transport have already been obtained. However, for the treatment of the proper spin 
current it is not sufficient to calculate f(k | s), it is rather necessary to treat V K /(fc,K | s) at K = 0. Therefore, we 
have to extend the perturbational approach with respect to k. For k^O and E = 0, we obtain 



f K0 = (k ■ k) 



Si- 



has 



(A7) 



- ih 2 n' in auJ K - 2ui k x lo k + AuJ k (uj k • cD K )/tr 

/ K o = (« • k)uj k 2 , 

m as s a + Alot. 



(A8) 



where the second term on the right-hand side of Eq. (A. 7) can be neglected. Finally, we compute the quantity that 
determines the spin-Hall current in Eq. 10, namely the contribution, which is proportional to k and E 
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aR z KE - 2(a) fc x R KE y 
a 2 + Auj 2 



with 
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R k e = — (k • k)f 0E - iuJ K {k)f 0E T-^T^/kO + —p. — htikfnE- 



(A9) 



(A10) 



After a lengthy but straightforward calculation, we obtain the result for the spin-Hall current in Eq. (|10|l from 
Eqs. (A.9), (A.10), and ©. 
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